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ABSTRACT 


A  method  utilizing  vector  representation  is  investigated  for  deter- 
mining a  faulty  element  in  passive  and  active  networks  by  simple 
external  measurements. 

A  large  system  may  be  considered  as  an  interconnection  of  a  number 
of  subnetworks.  By  utilizing  the  relationships  between  the  magnitudes 
of  a  transfer  function  at  various  frequencies  and  the  deviations  of  a 
circuit  element,  the  fault  simulation  curves  can  be  drawn.  The  fault 
identification  regions  are  defined  from  the  fault  simulation  curves. 
A  fault  identification  matrix  is  constructed  corresponding  to  the  de- 
fined fault  identification  regions. 

The  fault  identification  matrix,  when  premulti plied  by  a  vector 
whose  components  are  measured  from  a  network,  yields  another  vector 
whose  components  identify  a  network  element  which  is  faulty. 

A  test  procedure  for  the  fault  identification  method  is  presented 
and  verified  by  experiments. 
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I.   INTRODUCTION 

Recent  advances  in  the  design  of  electronic  systems  have  resulted  in 
an  ever-increasing  complexity  in  such  systems.  Despite  the  advances  in 
component  reliability  and  preventive  maintenance,  it  appeared  that  a 
great  many  unpredictable  or  random  failures  do  occur.  Such  complexity 
in  modern  equipments  used  by  the  industry  and  the  military  has  resulted 
in  a  serious  problem  of  providing  adequate  maintenance.  One  of  the  most 
difficult  tasks  is  the  location  of  malfunctions  which,  until  recently* 
has  been  left  entirely  in  the  hands  of  technicians. 

Some  techniques  for  fault  identification  have  been  studied  by  a 
number  of  investigators  [1,2,3,5,6,7,8,9,10]  and  a  computer  program  for 
fault  simulation  of  passive  linear  networks  without  mutual  inductance 
has  been  given  in  Ref.  [2]. 

In  this  paper,  the  problem  of  fault  identification,  which  requires 
the  least  number  of  external  test  points  of  the  system,  is  investigated 
and  a  procedure  which  will  automatically  identify  a  fault  or  faults  in 
passive  and  active  linear  networks  is  presented.  Toward  achieving  this 
ultimate  goal,  the  sensitivity  of  a  network  function  due  to  element  devi- 
ations is  "studied.  The  method,  which  requires  matrix  operations,  is 
introduced.  In  this  method,  the  fault  identification  regions  at  a  set 
of  selected  test  frequencies  are  defined  depending  on  the  network  sensi- 
tivity and  the  fault  identification  matrix  is  generated.  A  vector  is 
then  obtained  by  premulti plying  the  fault  identification  matrix  by 
another  vector,  and  the  resulting  vector  identifies  the  fault. 

The  single  element  fault  identification  method  is  discussed  and  the 
combined  element  fault  identification  method  is  introduced  by  utilizing 
the  single  element  fault  identification  method. 


The  computer  programs,  which  identify  a  fault  element  in  a  properly 
partitioned  network,  using  the  proposed  method  are  written  in  FORTRAN  IV 
and  included  in  Appendix  A.  The  implementation  of  the  computer  programs 
is  annexed  in  Appendix  B. 


II.   PRELIMINARY  CONSIDERATIONS 

A.   DEFINITION  OF  A  VECTOR 

Consider  a  two-port  network  with  the  external  terminals  and  variables 
as  defined  in  Fig.  2-1.  The  forward  voltage  ratio  transfer  function, 


Fig.  2-1. 
2-port  network. 

T12,  is  defined  as  the  ratio  of  the  output  voltage  to  the  input  voltage 
with  the  output  current  set  to  zero. 


T12  = 


I2=0 


(2-1) 


The  magnitude  of  T12,  designated  by  |T12|  ,  is  given  by  either  a 
function  of  the  frequency  with  fixed  element  values  or  a  function  of  an 
element  value  with  a  fixed  frequency.  That  is,  at  a  set  of  n  frequen- 
cies, with  fixed  element  values,  a  set  of  n  different  magnitudes  of  a 
transfer  function,  T12,  can  be  obtained.  Similarly,  corresponding  to  a 
set  of  k  deviations  of  an  element  value,  with  a  fixed  frequency,  a  set 
of  k  different  magnitudes  of  Tl 2  can  be  determined. 

Let  the  network  contain  only  the  following  types  of  passive  elements 
resistors,  inductors  and  capacitors.  Consider  an  element  which  has  a 
nominal  element  value  E  with  a  deviation  designated  by  aE.  ,  where  the 

J 

subscript  j  implies  the  j-th  deviation. 


With  a  change  in  element  value  aE.  of  E,  n  values  of  |T12|  can 
be  obtained  corresponding  to  n  different  frequencies.  These  values  of 
|T12|  are  designated  by 

IT12I.,   ,   |T12|.«  ,  —  ,   |T12|,,   ,  —  ,   1T12I. 

where  the  first  subscript  indicates  the  j-th  deviation  and  the  second 
subscript  indicates  the  i-th  frequency.  Next,  with  a  nominal  element 
value  of  E,  n  values  of  |T12|  can  be  obtained  corresponding  to  n  dif- 
ferent frequencies  and  designated  by 


T12|01   ,   |T12|02  ,  —  ,   |T12|01   ,  —  ,   |T12|0n 


The  normalized  value  of  |T12|..  is  defined  as 


iT12iji 


T12i0i 


where  i  =  1 ,  2,  —  ,  n. 

For  example,  with  a  deviation  aE.  of  E  and  a  set  of  5  frequencies,  the 

J 


5  values  of  the  normalized     |T12|,.     are 


IT12I  IT12I  IT12I 

I  '  '^hl  '       '12  i_ili 

m2 1  '  in?  I  »  '       IT12 1 


With  another  deviation,   aE     ,  of  E  and  the  same  set  of  5  frequencies, 
the  5  values  of  the  normalized     |T12L.     are 

|T12L.  |T12L«  |T12L, 


^2^0^      '         TT12T02      '  '      |T12|05 
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Let  a  vector  X  be  denoted  by 


X  =  (  x-j,  x2>  ---  ,  x.,  ---  ,  x  ,  c  )' 


(2-2) 


where  X  is  an  array  of  the  n  normalized  values  of  |T12|  with  respect 

to  a  set  of  n  frequencies  (with  a  deviation  of  aE.)  and  a  constant  c; 

|T12|..  J 

x,-  denotes  i-r-,0iJ1  (i  =  1,2,  ---  ,  n)  and  c  is  unity; 


T12 


Oi 


X  is  n  dimensional  because  the  last  component  is  always  unity; 

the  superscript  (')  means  a  transpose  of  an  array  of  X  or  a  matrix. 

B.  CONSTRUCTION  OF  FAULT  IDENTIFICATION  MATRIX  AND 
GENERATION  OF  VECTOR-Z 

Consider  a  simple  two-port  network  that  contains  a  resistor,  R,  and 

a  capacitor,  C,  as  shown  in  Fig.  2-2.  Let  the  nominal  value  of  R  be 

1  ohm  and  that  of  C  be  1  farad.  The  value  of  C  can  be  changed  from 


Fig-  2-2. 
2-port  RC  network. 

zero  to  infinity  and  3  nonzero  finite  values  of  C  are  indicated  as  shown 
in  Fig.  2-3. 


0.0 


0.5 


1.2 


1.5 


co 

— >- 


Fig.  2-3 


In  general,  consider  that  the  value  of  a  parameter  can  vary  from 
zero  to  infinity  and  this  range  of  values  can  be  divided  into  a  number 
of  subdivisions,  each  of  which  shall  be  called  the  j-th  deviation  of  the 
parameter  value,  where  j  =  1,  2,  —  .  For  example,  the  value  of  C 
(between  0  to  »)  is  divided  into  four  subdivisons,  that  is,  the  1-st 
deviation  is  0  <^C  <  0.5  ,  the  2-nd  deviation  is  0.5  ^  C  <  1.2  and 
similarly,  1.2  ^  C  <  1.5  and  1.5  ^C  represent,  respectively,  the 
3-rd  and  4-th  deviation.  • 

The  normalized  |T12|  is  given  by 


H!!ji  .  £U_  (,3) 


T12 


oi      yr+  (wc); 


The  values  of  the  normalized   ]T12|  corresponding  to  the  1-st  deviation 
through  the  4-th  deviation  are  shown  in  Fig.  2-4  according  to  a  set  of 
3  test  frequencies.  The  values  of  the  normalized  |T12|  corresponding 
to  the  j-th  deviation  where  j  =  1,2,3,4  are  best  seen  by  a  coordinate 
plane  whose  ordinate  represents  a  frequency  axis  and  abscissa  represents 
a  normalized  |T12|  axis  as  shown  in  Fig.  2-5. 

Let  a  "region-j",  designated  by  R-j ,  represent  an  area  on  the  fre- 
quency-normalized |T12|  coordinate  plane  which  corresponds  to  the  j-th 
deviation  of  an  element.  For  example,  R-l  represents  the  area  to  the 
right  of  line  A-0  (Fig.  2-5)  on  the  frequency-normalized  |T12|  coordi- 
nate plane  and  corresponds  to  the  1-st  deviation  of  C.  R-2  represents 
an  area  between  line  B-0  and  line  A-0,  including  line  A-0,  on  the  fre- 
quency-normalized |T12|  coordinate  plane  and  corresponds  to  the  2-nd 
deviation  of  C  as  shown  in  the  figure. 
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1-st 
deviation 


2-nd 
deviation 


0. 


w=l 


w=2 


w=4 


0.5 


1 > 


7  — * 


l3 — * 


1.265 


1.582 


1.845 


3-rd 
deviation 


4-th 
deviation 


1.2 


1.5 


0.906 

X 

0.861 

X 

> 

0.845 

X 

1  — > 


2  * 


3  > 


0.785 


-*1  + 


0.707 


0.68 


^3  > 


Fig.  2-4. 
The  normalized  |T12|  corresponding  to  the  j-th  deviation 
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Fig-   2-5. 
The  frequency-normalized    |T12|   coordinate  plane. 
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A  "normal  region"  is  defined  as  a  region  on  the  frequency-normalized 
|T12|  coordinate  plane  where  a  corresponding  element  deviation  is  in  a 
predefined  normal  range  of  deviation.  The  region  outside  of  this  normal 
region  is  defined  as  a  "faulty  region".  For  example,  if  the  2-nd  devi- 
ation is  a  predefined  normal  range  of  deviation,  R-2  is  a  normal  region. 
The  regions  outside  of  R-2,  that  is,  R-l ,  R-3  and  R-4  are  considered  to 
be  faulty  regions. 

Let  x.  be  a  variable  in  a  region  R-j  of  domain  x.  If  the  value  of 
x.  is  in  a  normal  region,  the  corresponding  element  is  considered  to  be 
normal.  For  example,  let  R-2  be  a  normal  region.  If  the  values  of  x. 
where  i  =  1,2,3  is  in  R-2,  the  element  C  is  considered  to  be  normal. 
Otherwise  the  element  with  the  values  represented  by  x.  is  considered 
to  be  faulty. 

Let  a  "boundary  point  value",  abbreviated  by  B.P.V.,  represent  a 
value  of  the  normalized  |T12|  which  lies  between  two  regions.  The 
boundary  point  values  of  a  region  may  or  may  not  belong  to  that  region. 
For  example,  the  boundary  point  values  between  R-l  and  R-2  lie  on  line 
A-0  and  these  values  do  not  belong  to  R-l  but  belong  to  R-2  as  shown 
in  Fig.  2-5.  The  3  boundary  point  values  corresponding  to  x-j ,  x2  and 
x3  in  R-l  are  1.265  ,  1.582  and  1.845  and  belong  to  R-2. 

A  "fault  identification  matrix",  designated  by  M,  is  defined  as  a 
matrix  whose  number  of  columns  represent  the  number  of  regions  and  the 
number  of  rows  represent  the  number  of  components  of  X  as  shown  in 
Fig.  2-6. 

A  vector  "  2   "  is  defined  as 

V     =  X'M  (2-4) 
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M  = 


(No. 

of 

regions  = 

No.  of  co 

lumns) 

R-l 

R-2  --- 

R-j  — 

R-r 

j_ 

._! 

i 

1 

"~T 

i_ 

._J. 

2 

-_T__ 

i 

n 

i_ 

._J 

n+1 

i 

i 

"7 

__"T 

(No.  of  X  components  =  No.  of  rov/.s) 

Fig.  2-6. 
Fault  identification  matrix. 


A  method  for  the  construction  of  M  and  the  generation  of  1   in  the 
n-dimensional  domain  x  with  r-regions  is  presented  as  follows.  The 
r-regions  R-l,  R-2,  ,  R-j,  ,  R-r  and  p  boundary  point  values 


by       b2 

2-7. 


,   b     ,  assigned  by  positive  numbers,  are  defined  in  Fig 


Corresponding  to  the  boundary  point  values,   the  r-regions  can  be 
arranged  with  respect  to  the  vector  component  x.   as  shown  in  Table  2-1. 
Let  the  regions,  arranged  as  shown  in  Table  2-1,  be  called  the  "identi- 
fication regions". 

A  fault  identification  matrix,  M,  will   be  constructed  corresponding 
to  the  defined  identification  regions.     For  every  B.P.V.   associated 
with  each  vector  component,  x.    ,  a  sub-fault  identification  matrix  is 
formed.     These  sub-fault  identification  matrices,  designated  by  m, ,  nu 
,  have  the  same  dimension  as  M. 
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R-l 


R-2 


R-j 


R-r 


k2  — 


xl- 


(P-D 


Fig.    2-7.       • 
n-dimensional   domain  x  with  r- regions 
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x2 
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bl 

<  x?  <  b? 

b,  <  x^  <  b? 

b?  <  x? 

1 

x. 
1 

1 

xn 

x  <  b 

n    p 

x  <  b 
n    p 

b  <  x„ 

P    n 

Table  2-1. 


Identification  regions 
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The  following  rules  are  given  for  the  construction  of  a  fault  identi- 
fication matrix. 
Rule  1: 

To  construct  a  sub-fault  identification  matrix  for  a  particular 
x.,  the  following  procedures  are  used. 

(a).      If  the  value  of  x.   in  R-j  is  greater  than  the  selected 
B.P.V.,   that  is  associated  with  x. ,  put  +1    in  the  i-th  row  of  j-th  column 
of  the  sub-fault  identification  matrix. 

(b).     If  the  value  of  x.   in  R-j  is  less,  than  the  selected 
B.P.V.,  that  is  associated  with  x. ,  put  -1   in  the  i-th  row  of  j-th  column 
of  the  sub-fault  identification  matrix. 

(c).     To  generate  the  entry  of  the  j-th  column  of  the  last 
row,  multiply  the  selected  B.P.V.   by  the  entry  value  of  the  j-th  column 
of  the  i-th  row  and  change  its  sign. 

(d).     All   entries  other  than  those  in  the  i-th  row  and  the 
last  row  are  set  to  zero. 
Rule  2: 

The  fault  identification  matrix,  M,   is  constructed  as  a  matrix 
sum  of  all- sub-fault  identification  matrices. 

Using  rules  1   and  2,  the  sub-fault  identification  matrices  and 
M  for  the  identification  regions  in  Table  2-1   are  constructed  as  follows. 
With  respect  to  x, ,  the  sub-fault  identification  matrices  are: 
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for  B.P.V.   of  b 


(P-1)' 


m. 


R-l  R-2 

-1.  -1. 

0.  0. 

0.  0. 

+b(p-D  +b(P-D 


R-J 


R-r 
+1. 

0. 

0. 


-b 


(P-1) 


1 
2 

i 

n 
n+1 


and  for  B.P.V.   of  b, 


m, 


R-l 

R-2 

+1. 

-1. 

0. 

0. 

0 
-b 


0. 


2  +b2 


R-J 


R-r 
+1. 

0. 

0. 
-bo 


i 

n 

n+1 


With  respect  to  x2,  the  sub-fault  identification  matrices  are 

—     R-j     — 


for  B.P.V. of  b   , 


m. 


R-l 

R-2 

0. 

0. 

-1. 

-1. 

0. 
+bo 


0. 
+bo 


R-r 

0. 

+1. 

0. 
-bo 


1 
2 
i 

n 
n+1 
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and  for  B.P.V.   of  b 


m, 


R-l  R-2 

0.  0. 

+1.  +1. 


R-j 


R-r 

0. 

+1. 


0. 
•b 


0. 


1  "1 


0 
-b 


1 
2 

i 
n 

n+1 


Finally,  with  respect  to  x   ,  the  sub-fault  identification  matrix  for 


B.P.V.   of  b     is 


in, 


R-l 

0. 
0. 


R-2  — -     R-j     -—  R-r 


0. 
0. 


-1.  -1. 


+b  +b 

P  P 


0. 
0. 

+1. 

-b_ 


1 
2 

i 

n 
n+1 


where  mf  means  the  last  sub-fault  identification  matrix, 
The  resulting  M  is  given  by 


M     =     m,   +  m«  +  itu  +  m,  + +  mf 


M     = 


R-l  R-2         —     R-j 

0.  -2. 

0.  0. 


-1.  -1. 


A1  A2 


R-r 
+2. 
+2. 

+1. 

A.. 


1 
2 

i 

n 
n+1 
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where  A,  =  b(p.1}  -  b]  +  —  +  bp  ,  A2  =  b(p.1}  +  2b2  -  b,  +  —  +  bp 

and  Ar  =  "  b(p-l)  "  2b2  "  bl  +  —  "  bp  • 

Next,  the  application  of  M  for  the  generation  of  Z  is  discussed. 

Let  the  normalized  values  of  |T12|  of  a  system  be  x, ,  x?, ,  x  . 

X  is  in  the  form  of 

X  =  J.  X-j  ,  Xp  >     »  x .  s  —  ,  x  ,  1.  J 
V    is  generated  as  follows. 

V       =    I'M  =  [  xr  x2, ,  xn,  1.  ]  M 

=  [  -xn+A]  ,  -2X]+ xn+A2  , ,  2Xl+2x2  --  +xn-:-Ar  ] 

J.   Z-.   ,         Zy      ,     ,   z»  ,         ,   z  J 


1   has  the  same  number  of  components  as  the  number  of  regions  and  z. 
serves  as  an  indicator  of  the  region  in  which  the  system  is  operating  as 
given  by  the  following  rule. 

Rule  3: 

If  a  system  is  operating  in  R-j ,  the  component  z.  of  Z  is  the  alge- 
braically  largest  component  in  Z. 

The  construction  of  M  and  the  generation  of  Z  are  illustrated  in 
the  following  examples. 

Example  2.1 .  Construct  M  and  generate  T  in  the  following  one 
dimensional  domain  x  with  three  regions. 

R-2  R-l  R-3 

5To  Toto 

x  <  5.0  5.0  <  x  <  10.0        10.0  <  x 
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Corresponding  to  the  given  boundary  point  values,  the  identification 
regions  are  as  follov/s. 

R-l  R-2  R-3 

x  5.0  £  x  <  10.0  x  <   5.0  10.0  <  x 

Using  rules  1   and  2,  the  sub-fault  identification  matrices  and  M  for  the 
identification  regions  are  constructed. 

The  sub-fault  identification  matrix  for  B.P.V.   of  5.0  is 

R-l  R-2  R-3 

+1.  -1.  +1. 

-5.0  +5.0  -5.0 


and  the  sub-fault  identification  matrix  for  B.P.V.   of  10.0  is 

R-l  R-2  R-3 

-1.  -1.  +1.  " 

+10.0  +10.0  -10.0 


The  resulting  M  is  given  by 


M    =    m,  +  m2 


R-l  R-2  R-3 

0.  -2.  2.  " 

5.0  15.0  -15.0 


Suppose  that  the  measured  value  of  x  is  13.0.  X  is  in  the  form  of 

t  -  r  x ,  i .  j 

[  13.0  1.  ] 
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Thus ,  V    is  obtained 


V     =  X'M  = 


[  13.0 

1.   ] 

0.            -2. 

2. 

_5.0         15.0 

-15.0 

E     5. 

-11. 

11.   ] 

[     z-,    , 

z2    . 

^ 

in  which  z3  is  the  algebraically  largest  component.  According  to  rule  3, 
the  system  is  operating  in  R-3.  The  measured  value,  of  x,  indeed,  belongs 
to  R-3. 

If  the  measured  value  of  x  is  6.0,  then 


V      =  X'M  =  [  6.0    1.  ] 


1.     -2.     2. 

5.0    15.0   -15.0 


[  5. 


3. 

z 


3.  ] 


2  •     Z3  ^ 


The  fact  that  z-,  =  5  ,  the  largest  of  all  the  entries  in  Z',  indicates 
that  the  system  is  operating  in  R-l.  ' 


20 


Example  2.2.  Construct  M  and  generate  Z  in  the  following  three 
dimensional  domain  x  with  four  regions. 


1.5 


2.0     2.5 


4.0 


5.0 


R-l 


R-2 


R-3 


R-4 


K  x, 


Xp  -=>! 


!<- 


<? —  X- 


■x,~ 


The  boundary  point  values  are  1.5  ,  2.0  ,  2.5  ,  4.0  and  5.0  .  The 
identification  regions  are  defined  as  follows. 


R-l 

R-2 

R-3 

R-4 

xl 

2.0  <  x,  <  4.0 

x1  <  2.0 

4.0  <  x. 

2.0 

<  x  <  4.0 

x2 

1.5  <  x2  <  2.5 

1.5  <  x2  <  2.5 

2.5  <  x2 

1.5 

<  x2  <  2.5 

x3  <  5.0 

x3<  5.0 

X3< 

5.0 

5.0 

<  X3 
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With  respect  to  x, ,  the  sub-fault  identification  matrices  are: 
for  B.P.V.   of  4.0, 


in. 


and  for  B.P.V.   of  2.0, 


m, 


R-l 

R-2 

R-3 

R-4 

-1. 

-1. 

1. 

-1. 

1 

0. 

0. 

0. 

0. 

2 

0. 

0. 

0. 

0. 

3 

4. 

4. 

-4. 

4. 

4 

R-l 

R-2 

.     R-3 

R-4 

1. 

-1. 

1. 

1. 

1 

0. 

0. 

0. 

0. 

2 

0. 

0. 

0. 

0. 

3 

-2. 

2. 

-2. 

-2. 

4 

With  respect  to  x?,  the  sub-fault  identification  matrices  are 
for  B.P.V.   of  2.5, 


m. 


R-l 

R-2 

R-3 

■  R-4 

0. 

0. 

0. 

0. 

1 

-1. 

-1. 

1. 

-1. 

2 

0. 

0. 

0. 

0. 

3 

2.5 

2.5 

-2.5 

2.5_ 

4 
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and  for  B.P.V.   of  1.5, 


m, 


R-l 

R-2 

R-3 

R-4 

0. 

0. 

0. 

0. 

1 

1. 

1. 

1. 

1. 

2 

0. 

0. 

0. 

0. 

3 

-1.5 

-1.5 

-1.5 

-1.5 

4 

With  respect  to  x3,  the  sub-fault  identification  matrix  for  B.P.V.  of 
5.0  is 


Plr 


R-l 

R-2 

R-3 

R-4 

0. 

0. 

0. 

0. 

1 

0. 

0. 

0. 

0. 

2 

-1. 

-1. 

-1. 

1. 

3 

5. 

5. 

5. 

-5. 

4 

The  resulting  M  is  given  by 


M  =  m-,  +  nu  +  nu  +  nu  +  nv 


R-l 

•     R-2 

R-3 

R-4 

0. 

-2. 

2. 

o.~ 

1 

0. 

o. 

•    2. 

0. 

2 

1. 

-1. 

-1. 

1. 

3 

8. 

12. 

-5. 

-2. 

4 

If  the  measured  components  of  X  are  x-,  =  3.0  ,  x2  =  2.0  and  x^  =  3.0  , 
X'  becomes 

X'  =  [  3.   2.   3.   1.  ] 


V     =     X' 


[3.  2.  3.  1.  ]  M  =  [  5, 


1.  ] 


L  zi  >  z9  ,  z-j  ,  zA   j 
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The  component  z,  is  the  albebraically  largest  component.  According  to 
rule  3,  the  system  is  operating  in  R-l.  Similarly,  if  the  measured 
components  of  X  are  x,  =  3.0  ,  x2  =  2.0  and  x~  =  6.0  ,  Z  is  given  by 

V     =  X'M  =  E  3.  2.  6.  1.  J  M  =  [  2.   0.  -1.   4.  ] 

1  '  ^2  *  ^3  '  ^4 

The  component  z»   indicates  that  the  system  is  operating  in  R-4. 
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III.  THE  PROPOSED  METHOD 

A.  PURPOSE 

The  purpose  of  the  proposed  method  is  to  investigate  a  feasible  method 
to  determine  whether  fault  detection  and  isolation  can  be  accomplished  by 
measuring  only  the  magnitude  of  the  system  transfer  frunction.  In  this 
method,  the  least  number  of  test  frequencies  and  test  points,  such  as  the 
normal  input  and  output  terminals  of  the  system,  will  be  used. 


B.  BASIS  FOR  THE  PROPOSED  METHOD 

Consider  again  the  2-port  network  shown  in  Fig.  2-1.  The  open  circuit 
impedance  functions  are  represented  by  the  equations  of  the  voltages  as 
linear  functions  of  the  currents  in  matrix  form  as  follows. 


Zll   Z12 


Z21   Z22 


1 


I 


(3-1) 


The  network  of  the  type  being  considered  can  be  externally  character- 
ized by  a  transfer  function  in  terms  of  open  circuit  impedances.  It  is 
easily  found  that  the  forward  voltage  ratio  transfer  function,  T12,  and 
the  reverse  voltage  ratio  transfer  function,  T21 ,  (defined  as  the  ratio 
of  the  input  voltage  to  the  output  voltage  with  the  input  current  set  to 
zero)  are 


T12  = 


:11 


(3-2) 


T21  = 


12 
:22 


(3-3) 


25 


In  general,  a  transfer  function  can  be  expressed  as  a  rational  func- 
tion of  the  complex  frequency,  s,  and  a  network  element,  E,  as  follows. 


m 

K    n 
i=l 

(s  +  z.) 

n 

n 
J-i 

(s  +  p.) 

T(s,E)  =  ^— (3-4) 


where  zero,  z.  ,  and  pole,  p.  ,  and  K  are  functions  of  E  and  m  and  n  are 
the  number  of  zeroes  and  poles  respectively. 

Taking  the  logarithm  of  Eq.  (3-4),  the  following  equation  is  obtained 

m  n 

lnT(s,E)  =  In  K  +  z     ln(s  +  z.)  -  z     ln(s  +  p.)  (3-5) 

1=1       ]    j=l       J 

And  taking  the  derivative  of  Eq.    (3-5),  the  following  equation  is 
obtained. 


d(lnTCs.E))    -   ^p 

t     Av  m        1       dz.  n         1       dp. 

=    1    dA    +       z 3L  .       z     _ 1        dE     (3-6) 

K     dE  1=1   s+z.     dE  j=T  s+p,     dE 


A  "transfer  function  sensitivity"  is  represented  by  Eq.  (3-6).  The 

transfer  function  sensitivity  is  closely  related  to  pole  and  zero  sensi- 

dp .       dz. 
tivities  where  y ■'  and  "tf^"  are  defined  as  pole  sensitivity  and 

zero  sensitivity  respectively.  In  other  words,  the  variation  of  T(s,E) 

may  take  the  form  of  a  change  in  location  of  a  pole  and/or  a  zero  of 

T(s,E). 
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From  Eq.  (3-4)  and  Eq.  (3-6),  it  can  be  seen  that  a  change  in  any 
one  of  the  pole-zero  configurations  will  cause  a  change  in  one  or  possi- 
bly both  of  the  voltage  ratio  transfer  functions.  Since  T12(s,E)  and 
T21(s,E)  depend  on  the  value  of  E,  it  is  not  hard  to  consider  how  a 
change  of  an  internal  element  is  related  to  a  change  of  an  external 
voltage  ratio  transfer  function. 

This  is  the  basic  philosophy  of  the  proposed  method  for  fault 
identification. 

C.  PROPOSED  TEST  PROCEDURE 

A  fundamental  fault  testing  procedure  for  the  proposed  method  is 
given  as  follows. 

Step  (1).  Partition  the  system. 

Step  (2).  Generate  a  symbolic  system  function. 

Step  (3).  Select  a  set  of  test  frequencies. 

Step  (4).  Simulate  a  fault. 

Step  (5).  Determine  the  fault  identification  regions. 

Step  (6).  Construct  the  fault  identification  matrix. 

Step  (7) .  Generate  X. 

Step  (8) .  Generate  Z. 

Step  (9).  Identify  the  fault. 

(a) .  Single  element, 
(b).  Combined  element. 

Step  (2),  step  (3)  and  part  of  step  (4)  have  been  studied  and  a 
computer  program  of  these  steps  for  the  passive  linear  network  without 
mutual  inductance  was  presented  by  J.  D.  Courville  in  Ref.  [2].  The 
detailed  discussion  of  these  steps  will  be  omitted.  The  other  steps  are 
discussed  in  detail  in  the  next  chapter  illustrative  examples. 
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IV.  NETWORK  FAULT  IDENTIFICATION 

A.  NETWORK  CONFIGURATIONS  AND  TRANSFER  FUNCTIONS 

As  illustrated  in  the  foregoing  discussions,  there  is  a  certain 
relationship  between  the  transfer  function  and  a  network  element  devia- 
tion. 

To  illustrate  the  proposed  method,  two  examples  for  passive  network 
configurations  are  presented. 

Example  4.1 : 


1/Eal  = 

Rl 

=  2  ohms 

1/Ea2  = 

L 

=  2  henry 

Ea3  = 

C 

'=  3  farads 

1/Ea4  = 

R9 

=  4  ohms 

Fig.  4-l(a) 
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Example  4.2: 


1/Ebl 

1/Eb2 

Eb3 

1/Eb4 


=  L1  = 

=  R  = 

=  C  = 

=  u  = 


1  henry 

1  ohm 

1  farad 

1  henry 


Fig.  4-l(b). 


For  example  4.1,  the  transfer  function  is 


T12     = 


F^LCs     +  (R1R2C+L)s  +  R-,   +  R2 


3s     +  6.5s  +  1.5 


K 


(s+0.262)(s+1.90) 


(4-la) 


where  K-,   is  1/3 
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For  example  4.2,  the  transfer  function  is 


L9s 

T12     =     £- 


L1L2Cs3  +  RL^s2  +   (Lj+L^s  +  R 


3         2 
s*5  +  s     +  2s  +  1 


K2s 


(s+0.57)(s2+0.43s+1.7623) 


(4-lb) 


where  K2  is  1    . 


The  network  semi  logarithmic  plots  for  the  magnitude  characteristic 
of  the  transfer  function  are  shown  in  Fig.   4-2a  for  Eq.  4-la  and  Fig. 
4-2b   for  Eq.    4-lb. 

B.      NECESSARY  DATA 

All   necessary  data  for  the  proposed  method  are  based  on  a  set  of 
chosen  test  frequencies. 

A  "fault  simulation  curve"   is  defined  as  a  response' curve  of  the 
normalized     |T12|     variations  due  to  the  simulated  deviations  of  an 
element  with  respect  to  a  chosen  test  frequency.     The  fault  simulated 
data  for     |T12|     can  be  obtained  from  the  computer  program  in  Ref.   [2]. 
The  fault  simulated  data  for     |T12|     and  the  normalized     |T12|     are 
tabulated  in  Tables  4-la  and  4-lb  for  Eq.   4-la  and  Eq.  4-lb  respectively. 

The  fault  simulation  curves  are  drawn  in  Graph  4-la-l   thru  Graph 
4-la-3  for  example  4.1    and  Graph  4-lb-l   thru  Graph  4-lb-3  for  example  4.2 
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C.  DETERMINATION  OF  FAULT  IDENTIFICATION  REGIONS 

The  proper  determination  of  the  identification  regions  is  the  most 
laborious  part  of  the  proposed  method. 

For  example  4.1,  a  set  of  3  test  frequencies  is  chosen  as  follows. 

w,  =  2.00  radians/sec 
Wp  =  1 .02  radians/sec 
Wo  =  0.02    radians/sec 

For  example  4.2,  a  set  of  3  test  frequencies  is  chosen  as  follows 

W-,  =  3.20  radians/sec 
w2  =  0.70  radians/sec 
w3  =  0.05    radians/sec 

For  a  test  frequency,  there  are  eight  basic  identification  regions, 
four  normal  regions  and  four  faulty  regions,  because  four  elements  are 
contained  in  the  examples.  That  is,  each  element  has  two  regions,  a 
normal  region  and  a  faulty  region.  The  four  normal  regions  are  not  the 
concern  of  the  propsed  method.  The  four  faulty  regions  may  be  sub- 
divided into  many  faulty  regions  to  identify  a  fault.  These  subdivisions 
of  the  faulty  regions  depend  on  the  transfer  function  sensitivity.  The 
subdivided  faulty  regions  are  called  the  fault  identification  regions. 

The  normal  regions  for  the  networks  are  defined  as  in  Table  4-2a  for 
example  4.1  and  Table  4-2b  for  example  4.2  according  to  the  fault  simu- 
lation curves.  The  regions  outside  of  the  normal  regions  are  considered 
to  be  faulty  regions.  The  fault  identification  regions  are  defined  as 
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Eal 


Ea2 


DEV 

MAG-T1 

MAG-T2 

MAG-T3 

N0R-T1 

N0R-T2 

N0R-T3 

(percent)  (w=2.0) 

(w=1.02) 

(w=0.02) 

1    90 

0. 11247 

0.26665 

0.79074 

1.8795 

1.8199 

1.1896 

70 

C.1CC92 

0.24123 

0.77166 

1.6865 

1.6464 

1.1609 

50 

0.C8928 

0.215C4 

0.74875 

1.4920 

1.4677 

1.1264 

30 

0. C7756 

0.18612 

0.72C74 

1.2961 

1.2840 

1.0843 

10 

0. C6576 

0.  16054 

0.68572 

1.0990 

1.0957 

1.0  316 

-10 

0.C5390 

0.13235 

0.64C69 

0.9008 

0.9033 

0.9639 

-30 

0.04199 

0.10262 

0.58067 

0.7017 

0.7072 

0.8736 

-50 

C. 03003 

0.07442 

0.49673 

0.5018 

0.5079 

0.7473 

-70 

0.C1603 

0.04485 

0.37120 

0.3014 

0.3  061 

0.5584 

-90 

C. 00602 

0.015C0 

0.16370 

0.1005 

0.1024 

0.2463 

2    90 

0.07451 

0.15650 

0.66464 

1.2452 

1.0681 

0.9999 

70 

0.07267 

0.15555 

0.66465 

1 .  2 1 44 

1.0617 

0.9999 

50 

C.C7C24 

0.15418 

0.66466 

1.1737 

1.0523 

0.9999 

30 

0.C6699 

0.15211 

0.66467 

1.1195 

1.0382 

0.9999 

10 

0.C626C 

0.14885 

0.66470 

1.0461 

1.0160 

1.0000 

-10 

0.C5662 

0.14347 

0.66473 

0.9462 

0.9792 

1.0000 

-30 

0.04851 

0.134C4 

0.66477 

0.8106 

0.9149 

1.0001 

-50 

0.03774 

0.11667 

0.664  84 

0.6306 

0.7963 

1.0002 

-70 

0.C2413 

0.C6429 

0.66497 

0.4032 

0.5760 

1.0004 

-90 

0.CC830 

0.C3149 

0.66482 

■  0.1386 

0.2149 

1.0002 

3    90 

0.03131 

0.C7739 

0.65942 

0.5232 

0.5282 

0.9920 

70 

O.C3502 

0.08652 

0.66C88 

0. 5853 

0.5905 

0.9942 

50 

0.C3974 

C.C9SC5 

0.66217 

0.6  641 

0.6692 

0.9962 

30 

0.C4591 

0.113C8 

0.66331 

0.7673 

0.7718 

0.9979 

10 

0.C5435 

C. 13341 

0.66428 

0.9C82 

0.9106 

0.9994 

-10 

0. C6656 

0.16  238 

0.66509 

1.1122 

1.1082 

1.0006 

-30 

0.C8576 

0. 2C620 

0.66574 

1.4331 

1.4094 

1.0015 

-50 

C. 12017 

0.28008 

0.66621 

2.0081 

1.9116 

1.0023 

-70 

0.19776 

0.41447 

0.66651 

3.3  047 

2.8288 

1.0027 

-90 

C. 45038 

0.61200 

0.66  66  5 

7.5262 

4.1770 

1.0029 

4    90 

0.C3991 

0.123C1 

0. 51212 

0.6669 

0.8395 

0.7704 

70 

O.C^-336 

0.12829 

0.53966 

0.7246 

0.8756 

0.8119 

50 

0.C4733 

0.13364 

0.57032 

'  0.7910 

0.9121 

0.8560 

30 

0.C5188 

0.  12694 

0.60467 

0.8670 

0.9433 

0.9097 

10 

0.C5704 

0.14407 

0.64342 

0.9532 

0.9833 

0.9680 

-10 

0.C6277 

0.14886 

0.68746 

1.0489 

1.0160 

1.0342 

-30 

0.06884 

0.15313 

0.73795 

1.1504 

1.0451 

1.1102 

-50 

0.07478 

0.15668 

0.79642 

1.2497 

1.0694 

1.1981 

-70 

C.C7973 

0.15932 

0.86491 

1.3324 

1.0874 

1.3012 

-90 

0.C8264 

0.16091 

0.94624 

1.3809 

1.0982 

1.4235 

Table  4-la. 

The  magnitudes  and  the  normalized  magnitudes  of 

T12  due  to  each  element  deviations  for  Example  4.1. 

(  The  nominal   magnitudes  are  0.05984  at  w=2.0   , 
0.1465  at  w=1.02  and  0.6647  at  w=0.02  ) 
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Ebl 


Eb2 


Eb3 


Eb4 


DEV 

MAG-T1 

fAG-T2 

MAG-T3 

NCR-T1 

N0R-T2 

N0R-T3 

(percent)  (w=3.2) 

(w=0.7) 

(w=0.05) 

90 

0.  24394 

0.6C490 

0.04991 

2.1299 

1.0142 

1.0009 

70 

0. 21254 

0.60589 

0.04  9  90 

1.8558 

1.0158 

1.0008 

50 

C. 18275 

0.6C623 

0.04990 

1.5957 

1.0164 

1.0007 

30 

0.15445 

0.60512 

0.04  9  89 

1 . 3^85 

1.0146 

1.0005 

10 

0.12751 

C.6C085 

0.04  9  67 

1.1134 

1  .0074 

1.0002 

-10 

0. 10166 

0.58945 

0.04965 

0. 6893 

0.9883 

0.9997 

-30 

C.C7739 

0.5615C 

0.C4960 

0.6757 

0.9414 

0.9987 

-50 

0.C54C3 

0.49499 

0.04  96  8 

0.4  717 

0.8290 

0.9963 

-70 

0.C3170 

0.35117 

0.C4  92  5 

0.2768 

0.5888 

0.9877 

-90 

0.01C34 

0.121C4 

0.C4468 

0.0903 

0.2029 

0.8961 

90 

0. 11935 

G. 64187 

0.09353 

1.C420 

3  .0762 

1.8759 

70 

0. 11886 

0.63708 

0.08398 

1.0378 

1.06  81 

1.6843 

50 

C. 11818 

0.63C43 

0.C7433 

1.0316 

1.0570 

1.4908 

30 

0.11718 

0 .62C85 

0. C6460 

1.0231 

1.040$ 

1.2956 

10 

0. 11563 

0.6C646 

0.05479 

1.0096 

1.0168 

1.0989 

-10 

0.11309 

0.56365 

0.04492 

0.9874 

0.9786 

0.9008 

-30 

0. 10852 

0.54525 

0.03499 

0.9475 

0.9142 

0.7018 

-50 

0.09938 

0.47653 

0.02502 

C. 8677 

0.7990 

0.5019 

-70 

C. C7892 

0.34967 

0.01503 

0. 6891 

0.5863 

0.3014 

-90 

0.C3330 

0.13439 

0.C0501 

C.2908 

0.2253 

0.1005 

90 

0.C5440 

0.93146 

0.C4997 

0.4749 

1.5617 

1.0022 

70 

0.06158 

0.63950 

0.C4995 

0.5377 

1.4075 

1.0017 

50 

C.07C95 

0.75  7  30 

0.C4992 

0.6195 

1.2697 

1.0012 

30 

C. 08369 

0.665  70 

0.C4990 

0.7307 

1  .1497 

1.0C07 

10 

0. 1C20C 

0.62398 

0.C4987 

0.8906 

1.0462 

1.0002 

-10 

0. 13C56 

0.57C89 

0.04985 

1.1400 

0.9572 

0.9998 

-30 

0.18130 

0.52510 

C.C4982 

1.5830 

0.8804 

0.9993 

-50 

0.29628 

0.46544 

0. C4980 

2.5869 

0.8139 

0.9988 

-70 

0.7985C 

C. 45090 

0.04677 

6. 9719 

0.7560 

0.9782 

-90 

1. C2460 

0.42067 

0.04975 

8.9460 

0.7053 

0.9978 

90 

0. 10913 

0.3596C 

0.02630 

0.9528 

0.6032 

0.5274 

70 

0.11115 

0.3  94  72 

0.C2  93  8 

0.9705 

0.6618 

0.5893 

50 

0.11275 

C.437C7 

C.C3  32  9  ' 

0. 9844 

0.7328 

0.6677 

30 

0.11386 

0.48950 

0.C3  840 

0.9942 

0.8207 

0.7701 

10 

0. 11445 

0.556C2 

0.04535 

0. 9993 

0.9322 

0.9095 

-10 

0. 11448 

0.643C7 

0.C5537 

0.9995 

1.0782 

1.110  5 

-30 

0.11395 

0.76155 

C.C71C6 

0.9949 

1.2768 

1.425  2 

-50  ■ 

0. 11289 

0.93135 

0.C9911 

0.9856 

1.5  615 

1.9  877 

-70 

0.U134 

1.19190 

0.16324 

0.9721 

1.9984 

3.2  740 

-90 

0.10935 

1.62772 

0.43908 

0.9548 

2.7291 

6.80  61 

Table  4-lb. 

The  magnitudes  and  the  normalized  magnitudes  of 

T12  due  to  each  element  deviations  for  Example  4.2. 

(  The  nominal   magnitudes  are  0.11453  at  w=3.2   , 
0.5964  at  w=0.7  and  0.04986  at  w=0.05  ) 
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in  Table  4-3a  for  example  4.1   and  Table  4-3b  for  example  4.2.      In  Tables 
4-3a  and  4-3b,  an  arrow  "+"  means  "more  than"  and  an  arrow  "4-"  means 
"less  than". 

D.      CONSTRUCTION  OF  FAULT  IDENTIFICATION  MATRIX 

The  sub-fault  identification  matrices  and  the  fault  identification 
matrix,  M,   for  the  fault  identification  regions  defined  in  Table  4-3a, 
are  constructed  step  by  step  according  to  the  rules  given  previously. 
The  boundary  point  values  for  x.  where  i   =  1,2,3  are  tabulated  in  Table 
4-4. 


B.P.V.'s 

xl 

1.2 

1.1 

0.9 

0.8 

x2 

1.2 

1.05 

0.95 

0.8 

x3 

1.06 

1.02 

0.98 

0.92 

Table  4-4. 


The  sub-fault  identification  matrices  with  respect  to  x-,   are  repre- 
sented by  m, ,  m?,  m^  and  m.    ,  the  sub-fault  identification  matrices  with 
respect  to  x?  are  represented  by  nv,  nv,  nv  and  nv  and  the  sub-fault 
identification  matrices  with  respect  to  x~  are  represented  by  riv,  m^, 
m,,   and  m,?.     All   of  the  sub-fault  identification  matrices  and  M  are 
given  as  follows. 
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E.  GENERATION  OF  X 

The  components  of  X  can  be  obtained  by  the  measurements  of  the 
normalized  |T12|  corresponding  to  the  chosen  test  frequencies. 

Suppose  that  it  is  desired  to  determine  which  element  is  faulty  in 
a  network  accessible  only  at  the  input  and  output  terminals.  The 
diagnosis  of  the  network  malfunction  can  be  accomplished  using  the  fault 
identification  matrix. 

The  value  of  X  will  fall  in  one  of  the  defined  fault  identification 
regions. 

F.  GENERATION  OF  Z  AND  FAULT  IDENTIFICATION 

The  generation  of  I   and  the  fault  identification  are  best  seen  by 
the  following  three  cases  for  example  4.1. 

Case  1.  Let  the  measured  values  of  the  components  of  X  at  the 
selected  test  frequencies  w, ,  w2  and  w3  be 

x,  =  .6 
x2  =  .7 
x3  =  .9 

1   is  generated  as  follows. 
V     =  X'M 


=  [  .6   .7   .9   1.  ]  M 
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z    = 


-2.4  -2.8  -3.6  +11.98 

2.4  +2.8  +3.6  +11.98 

-2.4  -1.4  +1.8  +  4.54 

-2.4  -2.8  +1.8  +  6.14 

1.2  +1.4  -1.8  -   1.06 

2.6  +1.4  -1.8  -   3.46 

-2.4  -2.8  -1.8  +10.14 

2.4  +2.8  +1.8  -  9.86 

-1.2  -1.4  -3.6  +  8.78 

-2.4  -1.4  -3.6  +10.38 

1.2  +1.4  +3.6  -   7.18 

2.4  +1.4  +3.6  -   9.58 


= 

3.18" 

-3.18 

2.54 

2.74 

-0.26 

-1.46 

3.14 

-2.86 

2.58 

2.98 

' -0.98 

_  -2.18 

10 


11 


12 


Case  2.     If  the  measured  values  of  the  components  of  X  at  the 
selected  test  frequencies  are 

x1   =  1.3 

x2  =  1.1 

x3  =  0.95 

V     =     X'M 

=     [1.3     1.1      .95     1.   ]  M 

=     [  -1.42         1.42         -.96       -1.56         1.84         2.04 

-1.36         1.64  .18         -.82         1.42         1.62  ] 


=     E       z. 


10 


z, 


z12  ] 
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Case  3.  If  the  measured  values  of  the  components  of  X  at  the 
selected  test  frequencies  are 


x1  =  2. 


x2  =  1.1 


x3  =  .98 


V     =  X'M 

=  [  2.    1.1  .98    1.  ]  M 

=  [  -4.34   4.34  -3.7    -4.3    3.18    4.78 

-4.22  4.5    -1.34  -3.74    2.94    4.54  ] 


C   Zl     z2     z3     z4     z5     z6 


Z7     Z8     z9     Z10     Zll     z12  ] 


From  the  defined  fault  identification  regions,  Table  4-3a,  the 
following  facts  are  known  for  each  case. 

Case  1.  The  algebraically  largest  component  in  I   is  the  component  z, , 
The  network-  is  operating  in  R-l  which  means  Eal  is  deviated  less  than 
-20  percent  from  the  nominal  value.  In  example  4.1,  R,  is  out  of  toler- 
ance. 

Case  2.  The  component  zg  in  1   is  the  algebraically  largest  component, 
The  netv/ork  is  operating  in  R-6.  The  element  Ea2  is  deviated  more  than 
60  percent  from  the  nominal  value.  In  example  4.1,  L  is  out  of  tolerance. 

Case  3.  The  component  zfi  in  T  is  the  algebraically  largest  component 
again  but  the  noticeable  fact  is  that  the  value  of  the  component  z,   is 
different  from  that  of  the  previous  case.  The  relative  value  of  the 
component  of  Z  may  suggest  the  degree  of  violation  of  the  faulty  element. 
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The  method  discussed  in  the  foregoing  illustrations  is  called  the 
single  element  fault  identification  method.  A  combined  element  fault 
identification,  abbreviated  by  C.  E.  F.  !•;,  method  is  a  method  in  which 
more  than  one  element  fault  is  identified. 

The  algorithm  of  the  single  element  fault  identification  method  is 
suitable  for  computer  implementation.  If  the  fault  identification 
regions  are  properly  defined,  then  the  construction  of  M  and  the  genera- 
tion of  Z  are  matrix  operations.  The  fault  identification  is  the 
comparison  of  the  values  of  the  components  in  Z. 

A  program  for  the  single  element  fault  identification,  designated  by 
program  1,  is  written  as  shown  in  Appendix  A  and  its  implementation  is 
annexed  in  Appendix  B.  The  results  of  the  computer  applications  for 
example  4.1  are.  shown  in  Fig.  4-3,  Fig.  4-4  and  Fig.  4-5.  The  same 
procedures  are  used  for  example  4.2  and  the  results  of  the  computer 
applications  are  shown  in  Fig.  4-6,  Fig.  4-7  and  Fig.  4-8. 
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DEFINED  BCUNDARY  PCINT  VALUES 
NUMBER  CF  FREQUENCIES  NV=  3 
NUMBER  CF  REGIONS       NR=  12 

LOWER  B.P.V.    UPPER  B.P.V. 


0.0 
0.0 
0.0 

0.80000 
0.80000 
0.92000 

1.20000 
1.20000 
1.06000 

100.00000 
100.00000 
100.00000 

o.c 

0.80000 
1.02000 

0.80000 
0.95000 
1.06000 

0.0 
0.0 
1.C2000 

0.80000 
0.80000 
1.06000 

1.10000 
1.C5000 
0.92000 

1.20000 
1.20000 
0.98000 

1.20000 
1.05000 
0.92000 

100.00000 
1.20000 
0.98000 

0.0 
0.0 
0.92000 

0.80000 
0.80000 
0.98000 

1.20000 
1.20000 
1.02000 

ICO. 00000 

100.00000 

1.06000 

o.eoooo 

0.80000 
0.0 

'  0.90000 
0.95000 
0.92000 

0.0 

0.80000 

0.0 

0.80000 
0.95000 
0.92000 

1.10000 
1.05000 
1.C6000 

1.20000 

1.20000 

100.00000 

1.20000 
1.C5000 
1.06000 

100.00000 

1.20000 

100.00000 

Fig. 

4-3. 

Fault  identification  regions 
printed  by  computer  for  Example  4.1 
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DEFINED  BOUNDARY  POINT  VALUES 
NUMBER  CF  FREQUENCIES  NV=  3 
NUMBER  OF  REGIONS       NR=  13 

LOWER  B.P.V.    UPPER  B.P.V. 


0.0 

0.93000 

0.77000 

0.70000 
0.95000 
0.99000 

0.0 

0.62000 
0.77000 

0.70000 
0.93000 
0.99000 

0.0 
0.0 
0.77000 

0.70000 
0.82000 
0.99000 

1.35000 
1.02000 
1.01000 

100.00000 
1.05000 
1.30000 

0.70000 

0.82000 
0.0 

0.95000 
0.93000 
0.77000 

0.70000 

0.0 

0.0 

0.95000 
0.32000 
0.77000 

0.0 
0.0 
0.0 

0.70000 
0.82000 
0.77000 

1.C3000 
1.05000 
1.30000 

1.35000 

1.30000 

100.00000 

1.35000 
0.82000 
0.77000 

100.00000 
0.93000 
0.99000 

1.35000 

0.0 

0.77000 

100.00000 
0.82000 
0.99000 

0.0 

1.20000 

1.01000 

0.70000 

100.00000 

1.30000 

0.95000 

1.20000 
1.30000 

0.98000 
100.00000 
100.00000 

0.95000 

0.0 

0.0 

0.98000 
0.82000 
0.77000 

Fig.  4-6. 
Fault  identification  regions 
printed  by  computer  for  Example  4.2, 
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V.      COMBINED  ELEMENT   FAULT   IDENTIFICATION 

Multi-element  fault  identification,  abbreviated  by  M.E.F.I.,  is  a 
somewhat  complicated  task.     However,   if  multi-element  fault  simulation, 
abbreviated  by  M.E.F.S.,   for  a  certain  system  can  be  made,  M. E.F.I,   is 
possible. 

Consider  the  number  of  times  of  M.E.F.S.   necessary  for  a  M.E.F.I.. 
For  a  3-element  system,  using  the  abbreviation  for  the  combination  of  m 
elements  out  of  n  elements,  nCm,  the  number  of  times  to  be  simulated  is 
as  follows. 

3C1     =     3  —  one  element  fault  combinations 
3C2     =     3  —  two  elements  fault  combinations 
3C3     =     1   —  three  elements  fault  combination 

Including  ewery  element  deviation  of  10  times,  each  element  deviation 
from  -90  percent  to  +90  percent  of  the  nominal   value  with  20  percent  of 
step  size,  the  number  of  fault  simulations  necessary  for  a  M. E.F.I,    is 
as  many  as 

3C1     10  =         30  for  one  element  fault  simulations 

3C2     10     10  =       300  for  two  elements  fault  simulations 

3C3     10     10     10     =     1000  for  three  elements  fault  simulations 

1330  for  total   number  of  fault  simulations 

Thus,  a  total  of  1330  multi-element  fault  simulations  are  required  even 
in  a  simple  3-element  system. 

The  number  of  possible  M.E.F.S.   required  for  a  M. E.F.I,   is  shown  in 
Table  5-1   and  it  becomes  astronomical   as  the  number  of  elements  increases 
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No. 

of  element 
2 
3 
4 
5 
6 

No. 

of  M.E.F.S 

110 

1330 

14640 

161050 

1771560 

Table 

5-1 

This  astronomical  increase  of  M.E.F.S.  discourages,  one  to  make  any 
further  investigation.  However,  another  method  can  take  its  place  as  is 
explained  as  follows.  A  complex  system  can  be  partitioned  properly  such 
that  the  probability  of  failure  involving  more  than  one  element  at  a 
time  is  usually,  small . 

In  this  section,  the  combined  element  fault  identification  (C. E.F.I.) 
method  is  introduced.  A  complex  network  shown  in  Fig.  5-1  is  partitioned 
properly  into  two  subnetworks  shown  in  Fig.  5-2.  Two  subnetworks  NW-1 
and  NW-2  are  regarded  as  independent  and  the  single  element  fault  identi- 
fication method  can  be  applied  to  each  network. 


Fig.  5-1. 
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Fig.  5-2. 

Suppose  that  the  fault  identification  matrices  for  NW-1  and  NW-2  are 
constructed  and  designated  by 

Ml     for   NW-1 
M2     for   NW-2 

Let  M  be  defined  as 


M 


Ml    '    0 

i 


0   .    M2 


(5-1) 


The  measured  values  of  the  components  of  X  can  be  obtained  with 
respect  to  a  set  of  selected  test  frequencies  for  NW-1  and  NW-2  independ- 
ently and  designated  as  follows. 


xir  xi2'  xi3' 


91  '    ??  '    /"\  * 


'  xln' 


'  x2m' 


for    NW-1 
for    NW-2 


where  n  and  m  are  the  number  of  selected  test  frequencies  for  NW-1  and 
NW-2  respectively.  Then  Z  are  generated  as  follows. 
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V     =     X'M     =     XT,  X2 


mi  ; 
o"; 


0 
M2 


=     Xl'Ml,  X2'M2 


=     t  zl     z2     z3  "—   zrd     z(rd+l)   ""  zr  ] 


where 


V 

II ' 

X2* 
IV 

Z2' 


XT  ,  X2' 

^  xll     x12     x13  xln  }t   1 

[  x21     x22     x23  x2m  1.   ] 

XTM1     =     [z,     z2  z3    -  zrd] 

X2'M2     =     [z(rd+1)  ------  zr     ] 


r  is  the  total   number  of  fault  identification  regions  for  both  networks 
and  rd  and  r-rd  are  the  number  of  fault  identification  regions  for  NW-1 
and  NW-2  respectively.     Z  will   indicate  two  faulty  elements.     One  from 
Zl   indicates  a  fault  in  NW-1   and  the  other  one  from  Z2  indicates  a  fault 
in  NW-2. 

This  is  the  C. E.F.I,   method  for  two  elements  fault  .identification. 
The  C. E.F.I,   method  can  be  extended  to  isolate  multi-element  faults. 

For  example,  suppose  that  the  network  shown  in  Fig.   4-1  a  is  replaced 
by  NW-1   and  the  network  shown  in  Fig.   4-1 b  is  replaced  by  NW-2  for 
C.E.F.I.   and  the  fault  identification  matrix  for  NW-1    is  Ml   and  the 
fault  identification  matrix  for  NW-2  is  M2.     The  resulting  M  is  constructed 
by  Eq.    (5-1)   and  shown  in  Fig.    5-3. 

Let  the  measured  values  of  the  components  of  X  be 
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XT  =  [  xn  x12  x]3  1.  ]  =  [  0.6   0.7  0.9   1.  ] 

X2'  =  [  x21  x22  x23  1.  ]  =  [0.6   0.94  0.8   1.  ] 

X'   =  XT,  X 2 '  =  [0.6   0.7   0.9   1.  0.6   0.94   0.8   1.] 

V   becomes 

T    =  X'M  =  XT  Ml,  X2'M2 

T    =  [  3.18   -3.18    2.54    2.74   -0.26  -1.46    3.14   -2.86 

2.58    2.98   -0.98  -2.18    3.89  3.87    3.63   -1.11 

3.61    3.37    3.57  -0.83   -0.15  -0.39    1.45   -0.45 
2.67  ] 


[  z]  z2  z3       z12  z13       z25  ] 


One  fault  identification  region  from  R-l  R-12,  represented  by 

z  z   respectively,  indicates  a  fault  in  NW-1  and  the  other  one 

fault  identification  region  from  R-13 R-25,  represented  by  z,3  -- 

—  z25  respectively,  indicates  a  fault  in  NW-2.  In  this  example,  the 

albegraically  largest  component  among  z-,  z,2  is  z-,  and  the  largest 

component  among  z,3 z25  is  z-,3.  The  fault  identification  region 

R-l  indicates  that  Eal  in  NW-1  is  deviated  less  than  -20  percent  from 
the  nominal  value.  The  fault  identification  region  R-13  represents  R-l 
in  NW-2  because  NW-1  has  12  fault  identification  regions  and  indicates 
that  Ebl  is  deviated  less  than  -28  percent  from  the  nominal  value.  Two 
more  examples  are  presented  in  Fig.  5-4. 

A  program  for  C.E.F.I.,  designated  by  program  2,  is  written  as  shown 
in  Appendix  A  and  its  implementation  is  annexed  in  Appendix  B. 
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VI.  EXPERIMENTS 


A.   PASSIVE  LINEAR  NETWORK  EXPERIMENT 

The  network  shown  in  Fig.  6-1  (identical  to  that  of  Fig.  4-1  a)  was 
used  in  the  experiment. 


I'o 


"1 


L 


-o 


<H 


-OZ 


Rl 

= 

2  ohms 

L 

= 

2  henry 

C 

= 

3  farads 

R2 

= 

4  ohms 

Fig.  6-1. 

Passive  network 

The  following  relationships  are  used  for  frequency  and  magnitude 
scaling: 


R* 

= 

b  R 

L* 

= 

b/a  L 

C* 

= 

1/ab  C 

f* 

= 

a  f 

where  "a"  is  the  frequency  scale  factor  (dimensionless ,  real  and 
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positive  constant),  "b"  is  the  magnitude  scale  factor  (dimensionless, 
real  and  positive  constant)  and  a  scaled  nominal  element  value  is 
identified  by  an  asterisk  superscript. 

With  a  =  10000  and  b  =  100  ,  the  scaled  nominal  element  values 
were 

1/Eal  =  R,  =  200  ohms 

1/Ea2  =  L  =  20  mill i  henry 
Ea3  =  C  =   3  micro  farads 
1/Ea4  =  R2  =  400  ohms 

and  a  set  of  3  selected  test  frequencies  became 

w-j  =  20000  radians/sec,  f]  =  3185  Hz 
w2  =  10200  radians/sec,  f^  =  1624  Hz 
w3  =    200  radians/sec,   f3  =  31.8  Hz 

The  experimental  results  for  fault  simulations  at  a  set  of  selected 
frequencies  were  tabulated  in  Table  6-1  for  the  network  shown  in  Fig. 
6-1.  Comparing  Table  6-1  with  Table  4-1  a  and  taking  into  consideration 
the  accuracy  of  the  measuring  devices,  it  is  concluded  that  the  experi- 
mental results  in  Table  6-1  and  the  theoretical  results  in  Table  4-1  a 
are  close. 

Consider  the  following  four  cases. 

Case  1.  If  the  resistor,  R, ,  was  changed  to  130  ohms  from  the 
nominal  value  of  200  ohms,  the  magnitudes  of  T12  were  0.0935  ,  0.222 
and  0.755  at  a  set  of  3  selected  test  frequencies  f-,  ,  f^  and  f^ 
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respectively.     Since  the  nominal  magnitudes  of  T12  were  the  same  as  used 
in  Section  IV-C,   the  normalized  magnitudes  of  T12  became 


x1     =     1.562 


x2     =     1.515 


x3     =     1.138 


Case  2.     If  the  inductor,  L,  was  changed  to  30  milli  henry  from  the 
nominal   value  of  20  milli   henry,  the  magnitudes  of  T12  at  a  set  of  3 
selected  frequencies  were     0.0462   ,  0.131   and  0.664.     The  normalized 
magnitudes  of  T12  were 


xl 

s 

0 

777 

x2 

= 

0 

.891 

x3 

= 

1 

.0 

Case  3.     If  the  capacitor,  C,  was  changed  to  4  micro  farads,  the 
magnitudes  of  T12  were     0.448   ,  0.109     and     0.664     and  the  normalized 
magnitudes  of  T12  were 


x]     =     0.748 


xz     =     0.743 


*3     =     0.98 


Case  4.  If  it  was  made  that  the  resistor,  R-,  changed  to  250  ohms, 
the  magnitudes  of  T12  were  0.0453  ,  0.13  and  0.56  and  the  normalized 
magnitudes  of  T12  were 
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X1  =  0.757 


x2  =  0 . 889 


x3  =  0.84 


Using  the  fault  identification  matrix,  Fig.  4-4,  Z  and  the  fault 
identification  were  generated  as  shwon  in  Fig.  6-1. 

From  Fig.  6-1  and  the  defined  fault  identification  regions,  Table  4-3a, 
the  following  facts  were  known  for  each  case. 

Case  1.  The  component  z2  was  the  largest  component  in  T  and  X 
belonged  to  the  fault  identification  region  R-2.  Eal  was  deviated  more 
than  20  percent  from  the  nominal  value.  R,  was  out  of  tolerance. 

Case  2.  The  network  was  operating  in  R-3.  Ea2  was  deviated  to 
somewhere  between  -30  percent  and  -50  percent  from  the  nominal  value. 
L  was  out  of  tolerance. 

Case  3.  The  network  was  operating  in  R-7  and  Ea3  was  deviated  more 
than  25  percent.  C  was  out  of  tolerance. 

Case  4.  Ea4  was  deviated  more  than  50  percent  which  was  indicated 
by  R-10. 

B.  ACTIVE  LINEAR  NETWORK  EXPERIMENT 

The  network  shown  in  Fig.  6-2  was  used  for  active  linear  network 
experiment.  The  frequency  response  of  the  active  network  was  drawn  in 
Fig.  6-3  using  the  computer  program  ECAP.  A  set  of  selected  test  fre- 
quencies was  chosen  as  follows. 

f-  =    10  Hz 

f2  =   100  Hz 

f3  =  1000  Hz 


70 


CO 
I 

a 


inaror-orovGr^-vOcM 
mm  cu  o  oust  c\j  <-o  ro 
0*ocMO^covor->j-inro 

t\iHHOO(Mcr-in(vj 


cr>o^o^a*o>>T<rvj->3-"J" 

cucocococoooooo 

^a>0'0>o>ooooo 


>r>j-<ro^o,.>r<J'vrcr 

H«fr-wcuajoooH 
ff-o>OHi'a>o>oooo 
o\a>^u^o^oooo 


cooinr-coinaoocor- 
>Tcocor^o^r-iu>r--r- 

h-  rH  vO  »-<  Is-  CO  rH  0~>  O  CM 

r-  co  co  a-"  cj^  OrH  r-i  ro  sr 


lOOOOO     ooooo. 


OOOOOOHHHH      OOOOOi 


eg 

I 

a: 
a 


I 

ex 
o 


CMO^COOCr>CMCMinvO 

oo  co  r-  o  cm  o  so  co  ro  r- 
cornvor^-o^oo^cro^o^ 

NvOg-NOO^vO^CvJO 


r~cj\rH>j-rHor-r— coo 
-sr-rHr-Ovor^-rHCMr- 
o  umt»  m  i-i  p-  o  o  in  o-> 

OOOOOO'O'MDH 


inooor-coo-iineor- 
inr-cMr~r^o^rorHcoo> 

c\jcoomo<-trr>mo-4" 
ir\ir\sor-Oi-»sj-o^r^rH 


r-  vor-r-  o  a-<N  voo  o> 
cm  co  r-  oo  o  o  vj  f-t  cm  co 

COCCO'0>0>00000 


HHHHHOOOOO       HHHHHOOOOO       OOOOO^HHtM<t       OOOOOt 


co  o  ro  ih  ct^  -4-  o  o  in  o 
>3-ocMu\r--cMrHror-rH 
o  cm  a-- coo  ooooo 

COvO<J-CMOO^>tfMO 


vOvOr-tfococMoooror*- 

r-iHM>OMf>fOirig-oo 

sJ-cMr^-rHin-g-rHrHCOrO 
CM  CM  rH  rH  O  O  i  CO  vQ  CO  rH 


OCMrHfOrHOOO^TrH 
co  oo  O  O  vi  •  co  co  cm  u  ^  ro 
rHf-OvOOCMsOCMCMO 

ininsor-OrHsj-omm 


rHvor^cocMinr-coo>o 

O  00  CO  CM  *4"  vj-  O  CX  rH  r- 

vOH<a<oinuMAir\«i-o 
vO  I*-  Is"  co  o  O  •-<  i  \)  co  co 


IOOOOO 


•  ooooo     ooooorHrHCMmr-    ooooo. 


o 

c 

<u 
> 


CO  ^-s  oooooooooo 

H-  "OOOOOOOOOO 

I  **  Oir\OinOr-iOmom 
O  +j  oi—mcMOvj  ooor-in 
<t    ^  cor-h-r-ovom-4-mrH 


oooooooooo 
oooooooooo 
^sj-sfvf-sj-ininininin 

vOvO^O^vO^OvOOnOvO 
vOvOO^vOnOvOnOvOvO 


oooooooooo 
oooooooooo 
OrHrOsf-^-  sfu"\ininvO 
ir,vOvO*o^oo-oovOvO 


oooooooooo 
oooooooooo 
mrH<foinoo^r-oo 
rH-<t-r— rHvj-ooroo-osi- 
inininvo^-OvOf-r^coo 


oooooooooo  oooooooooo  oooooooooo  oooooooooo 


<M  ^  oooooooooo 

I-      CMOOOOOOOOOrO 
I    **-  CMOinr-OcMCMrom*tf- 

O       ..     <£>  sf  rH  00  v£>  CO  O  r- vf  rH 


rrj 


oooooooooo 
oooooooooo 
vOin^-f\jco(<i(CivOr-ia> 
ininininsfs^cnrHcocM 

CMCMlMrHrHrHrHOOO         HHi-<HHHHr-iOO 


OOOOOOOOOO  oooooooooo 


oooooooooo 
c^  0  r- .h  n  i  >$- o  oo  o  co 

r~-COO>rHO>NOrHCOrHO 
OOOtHrHrHCMCMvfvO 


OOOOOOOOOO 

cMooroOcoo*cor~-o»H 
cm  cm  co  ro  >}- st- in  in  o  so 


oooooooooo  oooooooooo  oooooooooo  oooooooooo 


rH  ^^  oooooooooo 

h~     ,—ooroor-ooinaoin 

I    q-  ooh-o^vCin^  cMO^r-m 

<S>         oooor^-oin-4-cMrHO 

<     *f    rHOOOOOOOOO 
-«-     fo 


oooooooooo 
mrHCMcoocor-coom 
>l-rno^m-ooo»ococo 
r-r-h-vO^O'CNvS-mcMO 
oooooooooo 


oooooooooo 
OvommrHCMO'OOO 
rH<j  oinvj-r>-r>-rHtr»o* 
m  co  en  st-  u  \  ->o  oo  cm  u^  <r 

UOOOOOOHH«J- 


OOOOOOOOOO 
m  O  CM  O  rH  rH  si- 1  o  rf  1  o 

0>nihHhf<io>mor<"i 
co  <f  «4-  in  u\  vi)  vO  r~  co  co 
oooooooooo 


oooooooooo 

to 


oooooooooo  oooooooooo  oooooooooo 


E 
x: 
o 


CCJ         «>» 

<y   3    mr— co<*-cSCSv£>ov£>o 

s  ri  o  ^roinoQCMoocx^o 
«  r°  „~v-,_-<o(cxvfsoo 


lor-oO-^-fScS'^'     ^3 


o.— irr>i/->ooCSCo  o<sD  o 


h-  >-<  in  cy>  ro  r-  *h  in  cr>  m 
in  In  <t- en  rocS  cH  »h 


OL->sDOO^-Lo>HOCn;0 

cvj  rs  cS  cr%  rr>  <t  lnoo  r-o  o 


r — 

fO 
I- 


</) 

to 

fO 
Q_ 

s- 
o 
q- 

fD 

+-> 

fO 
"D 

•o 
a» 
+-> 


r3 


fl3 
+J 

C 
OJ 

E 
u. 

OJ 
O- 
X 


cc: 


-r-> 


>    Qj  oooooooooo 

LU       O    Or^inCOrHrHrOmh-O'' 

°    fe  I    I    I   I    l 

Q- 


oooooooooo 

o^r^incorHr-^comr^o 

I   I   I   I   I 


CO 


oooooooooo 

o  r~  in  c<^  rH  r-i  co  m  r~-  c 

I   I   I   I   I 


CM 


oooooooooo 

or-inrOrHrHcoir  r-a> 

I    l    I    I    I 


4-> 

C 

OJ  r— 

E  ro 

oj        uj 


CM 
LU 


CO 

ro 


ro 
UJ 


71 


Is- 


in 


o 


o* 

CM 


in 

ro 


vO 


o 

i- 

> 


00 


<\J 


m 

2 

C\J 

o 

t 

•— i 

o 

o 

I 

UJ 

CL 

z 

vO 

Q 

r- 

»—t 

• 

t- 

r- 

< 

O 

•— « 

LL 

H-« 

ro 

H- 

o 

Z 

• 

UJ 

r- 

o 

CO 


ro 


o 

•—I 

o 

UJ 


O     O 


o 


UJ 

O 


in 


o 

z: 

vO 

o 

• 

•— i 

*4" 

o 

U-i 

0l 

z 

ro 

o 

in 

•—1 

• 

h- 

r-* 

< 

o 

•— « 

UL 

t— 1 

CO 

H- 

•4- 

z 

• 

UJ 

nj 

D 

O* 
CO 


o 

■z. 

ro 

o 

• 

►— 1 

*fr 

CO 

UJ 

<X 

in 


(\j 


o 


ro    o 


4-> 

c 
oj 


i- 

a. 
x 

OJ 

s- 
o 


OJ 


2 

o 

1 — 

OJ 

> 

*— 1 

1 

to 

h- 

10 

00 

< 

ro 

o 

CJl 

Q. 

»— t 

•r— 

S- 

u_ 

U, 

o 

•— 1 

<+- 

I 

u 
o 

O 

OJ 


o 
o 


eg  r> 

ro  < 

•  a. 
o 


o 
o 


»-i    < 
•    u_ 


o 
o 


1-*    id 
•     u. 


o 
o 


r- 


ro 


IP 
u. 


o 


OJ 


ro 


in 


o>j- 

•    • 

cr>co 


in 


o 
o 


co 


ro 


vOao 
rOCT< 


co 


-4" 


Sj" 


COO'' 
ro<\J 

•  • 


-4" 

CO 


CO 


O 
O- 

• 
ro 


inh- 

•   • 


oc 

o 

o 

UJ 

> 


in 


>or^ 

•  • 
or- 
I 


co(\j 

O^CO 


co 


in 

• 
o 


h-cr> 
in^-t 

•  • 


(NJvO 

ino^ 


72 


H.  ^'ce©  (f) 


00 

K  ohms 

10 

K  ohms 

1 

K  ohms 

5 

K  ohms 

5 

yf 

50 

irf 

Ql  —  2N1414 


H.  =  1300  ohms 


1/H 


le 

oe 

1 
re 

fe 


36  K  ohms 
2.9  10~4 
42 


Cno  "  26   Pf 


oe 
GM 


Hre/Hie 


(b)  small  signal  ECAP 

equivalent  circuit  for  Ql 


Fig.  6-2. 
Active  network 

Since  there  v/as  no  computer  program  available  for  fault  simulations 
of  the  active  case,  the  necessary  data  were  obtained  from  the  experi- 
mental results.  The  fault  simulated  data  were  tabulated  in  Table  6-2. 
The  fault  simulation  curves  were  drawn  in  Graph  6-1,  Graph  6-2  and  Graph 
6-3.  The  normal  regions  were  defined  in  Table  6-3  and  the  fault  identi- 
fication regions  were  defined  in  Table  6-4.  The  fault  identification 
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matrix  was  generated  by  the  program  1    (in  Appendix  A)  and  is  shown  in 
Fig.   6-5. 

1  and  the  fault  identification  were  obtained  for  the  following 
three  cases. 

Case  1.     Let  the  resistor,  R2>  be  changed  to  14  K  ohms  from  the 
nominal   value  of  10  K  ohms.     The  measured  magnitudes  of  T12  were     11.71    , 
91.0     and     179.0     at  3  test  frequencies  f,    ,  f?  and  f~  respectively. 
Then  the  normalized  magnitudes  of  T12  became 


x1     =     1.03 


x2     =     1.08 


x3     =     1.26 


Case  2.  If  the  resistor,  R  ,  was  changed  to  7  K  ohms  from  the 
nominal  value  of  5  K  ohms,  the  measured  magnitudes  of  T12  were  15.9  , 
126.0  and  190.2  at  3  corresponding  test  frequencies.  The  normalized 
magnitudes  of  T12  were 

X]  =  1.39 

x2  =  1.5 

x3  =  1.34 

Case  3.  If  the  capacitor,  C2>  was  changed  to  35  micro  farads  from 
the  nominal  value  of  50  micro  farads,  the  measured  magnitudes  of  T12 
were  9.5  ,  75.0  and  143.0  at  3  corresponding  test  frequencies.  The 
normalized  magnitudes  of  T12  were 
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x,     =     0.834 


x2     =     0.892 


x3     =     1.007 


Using  M,  shown  in  Fig.   6-5,  1  and  the  fault  identification  generated 
by  the  computer  (program  1)  are  shwon  in  Fig.   6-6. 

From  the  defined  fault  identification  regions,  Table  6-4,  the  follow- 
ing facts  were  known  for  each  case. 

Case  1.     The  component  z3  was   the  algebraically  largest  component   in 

Z.     R-3  indicated  that  Rp  was  deviated  more  than  30  percent. 

Case  2.     The  fault  identification  region  R-8  indicated  that  R    was 

3  c 

deviated  more  than  20  percent. 

Case  3.  The  fault  identification  region  R-13  indicated  that  C?  was 
deviated  between  40  percent  and  70  percent. 

C.   PASSIVE  AND  ACTIVE  NETWORK  C. E.F.I.  EXPERIMENT 

Suppose  that  the  networks  shown  in  Fig.  6-1  and  Fig.  6-2  were 
considered  as  an  interconnection  of  two-port  network  consisting  of  two 
cascade  connected  subnetworks  and  regarding  Fig.  6-1  as  NW-1  and  Fig.  6-2 
as  NW-2,  The  fault  identification  matrix  for  C. E.F.I,  experiment  were 
shown  in  Fig.  6-7. 

Applying  the  same  procedures  as  in  the  C.E.F.I.  method  previously 
discussed,  the  network  malfunctions  could  be  isolated.  The  C.E.F.I.  was 
accomplished  for  three  cases  with  the  same  X  used  in  the  passive  and 
active  network  single  element  fault  identification  experiments.  The 
resulted  fault  identification  was  shown  in  Fig.  6-8. 
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DEFINED  BOUNDARY  POINT  VALUES 
NUMBER  OF  FREQUENCIES  NV=  3 
NUMBER  OF  REGIONS     .  NR=  14 

LOWER  B.P.V.    UPPER  B.P.V. 


0.80000 

0.0 

0.0 

0.92000 
0.80000 
0.66000 

0.0 
0.0 
0.0 

0.80000 
0.80000 
0.66000 

1.02000 
1.07000 
1.20000 

1.05000 

1.20000 

100.00000 

1.05000 
1.07000 
1.20000 

1.20000 

1.20000 

100.00000 

1.20000 
1.07000 
1.20000 

100.00000 

1.20000 

100.00000 

0.92000 
0.80000 
0.66000 

0.97000 
0.93000 
0.84000 

0.0 
0.0 
0.66000 

0.80000 
0.80000 
0.84000 

1.20000 
1.20000 
1.20000 

100.00000 
100.00000 
100.00000 

1.05000 
1.07000 
0.84000 

'1.20000 
1.20000 
0.99000 

1.20000 
1.20000 
0.84000 

100.00000 

100.00000 

0.99000 

1.20000 
1.20000 
0.66000 

100.00000 

100.00000 

0.84000 

0.92000 
0.93000 
1.01000 

0.97000 
0.97000 
1.20000 

0.80000 
0.80000 
1.01000 

0.92000 
0.93000 
1.20000 

0.0 

0.0 
1.01000 

0.80000 
0.80000 
1.20000 

Fig. 

6-4. 

Fault  identification  regions 
printed  by  computer  for  Active  network 
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VII.      FEASIBILITY  OF  FAULT  PREDICTION 

The  feasibility  study  of  fault  prediction  is  quite  complicated. 
However  a  brief  discussion  of  this  topic  is  included  in  this  section 
Suppose  that  M  of  a  certain  system  is 


M    = 


all     a12     a13 


a21     a22     a23 


anl     an2     an3 


Al       A2       A3 


dl3 


l2j 


nj 


llr 


l2r 


nr 


and  X'   is 


X'     =     [  x]     x2     x3     ---     xi     — -     *n     1.   ] 


T'   becomes 


V     =     X'M     =     [  x1     x2     xp     1.  ]  M 


■     E  z1     z2 


zj     -     zr] 


where 


Zj     =    Vlj  ♦  x2a2j  + 


+  x  a    .  +  A. 
n  nj         j 


j=l,2,—  ,r 


x.   is  positive  real   number  (  1=1,2,— ,n)  and  a,  .  and  A.  are  real 

I  J  J 

constants. 
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Since  z.  is  linear  functions  of  x.s,  the  variation  of  the  value  of 
z.  is  monotonic. 

If  the  system  is  faulty  and  z.  is  algebraically  the  largest  component, 
the  system  is  operating  in  R-j  and  an  element  which  is  out  of  tolerance 
can  be  isolated  from  the  predefined  fault  identification  regions. 

Taking  the  network  shown  in  Fig.  4-1  a  as  an  example  for  the  fault 
prediction,  there  are  8  boundary  point  values  between  the  normal  region 
and  the  faulty  region  at  a  test  frequency.  The  largest  components  in  Z" 
at  these  boundary  point  values  are  shown  in  Table  7-1.  From  this  table 


Element 
(percentage) 

Eal 

(-20) 

Eal 

(+20) 

Ea2 

(-30) 

Ea2 

(+25) 

Ea3 

(-18) 

Ea3 

(+25) 

Ea4 

(-20) 

Ea4 

(+25) 

1     A2    *3 

(at  boundary) 


(0.80 
(1.20 
(0.80 
(1.10 
(1.20 
(0.80 
(1.10 
(0.90 


0.80 
1.20 
0.95 
1.05 
1.20 
0.80 
1.05 
0.95 


0.91) 
1.06) 
1.02) 
0.98) 
1.02) 
0.98) 
1.06) 
0.92) 


TheJIargest  component 
in  Z  at  boundary 


1.90 
1.86 
1.48 
1.28 
1.78 
1.78 
1.36 
1.40 


Table  7-1. 

the  following  observations  are  made: 

1)  If  the  deviation  of  Eal  element  is  to  be  decreased  less  than 
-20  percent,  the  largest  component  in  Z  must  be  algebraically  greater 
than  1.90;  if  the  deviation  of  Eal  element  increased  more  than  +20 
percent,  the  largest  component  in  Z  must  be  greater  than  1.86.  The  same 
idea  can  be  applied  to  all  elements. 
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2)  If  z.  for  all  j  is  smaller  than  1.28,  the  smallest  value  in 
Table  7-1,  the  network  is  working  in  normal  condition. 

3)  If  z.  for  all  j  is  larger  than  1.90,  the  largest  value  in  Table 
7-1,  the  network  is  working  in  one  of  the  faulty  regions. 

4)  If  the  largest  component  z.  in  Z  is  in  between  1.28  and  1.90, 
the  network  may  or  may  not  be  in  a  normal  region.  For  instance,  if  the 
largest  component  z.  is  1.50  and  the  network  is  still  in  normal  condition; 
then  it  is  known  that  Ea2  and  Ea4  are  not  faulty  and  Eal  or  Ea3  may  begin 
to  approach  one  of  the  faulty  regions.  If  the  largest  component  z.  is 
1.50  and  the  network  is  in  faulty  condition,  then  it  is  known  that  Ea2 

is  faulty. 
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VIII.  CONCLUSION 

A  fault  identification  technique  has  been  investigated  in  this  paper. 
Although  the  examples  used  are  simple,  the  fault  identification  proce- 
dures presented  here  may  be  extended  to  more  complex  systems. 

The  theoretical  basis  of  the  proposed  method  is  described  and  the 
technique  has  been  applied  to  both  passive  and  active  networks. 

The  main  differences  in  the  approach  used  in  this. paper  for  fault 
identification  over  the  gain  signature  approach  cited  in  References  [3], 
[8]  and  [10]  are  summarized  as  follows. 

1).  Fault  identification  regions  are  defined  uniquely  according  to 
the  fault  simulation  curves  which  are  constructed  from  the  system  sensi- 
tivity instead  of  assignment  of  gain  signature. 

2).  Direct  application  of  the  measured  values  for  fault  isolation 
can  be  accomplished  instead  of  gain  signature  conversion. 

3).  Vector-Z  indicates  an  exact  faulty  element  automatically  rather 
than  a  gain  signature  comparison. 

4).  The  selection  of  a  set  of  test  frequencies  is  not  critical  as 
in  the  gain  signature  approach. 

5).  By  applying  the  matrix  operation  to  the  fault  identification 
matrix  of  the  single  element  fault  identification  method,  one  can  detect 
not  only  a  fault  of  the  one-element  type  but  also  a  multi-element  fault 
in  the  given  system. 

The  proposed  method  uses  only  one  network  function.  However,  it 
gives  no  ambiguities  as  demonstrated  in  the  examples. 
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APPENDIX  B 
IMPLEMENTATION  OF  THE  COMPUTER  PROGRAMS 


INPUT  DATA 


PROGRAM  1 

Card 
number 


Card 

column 


Information 


Form 


NV+1 


NV+2 


--  10 

-  20 

--  10 

--  20 

--  10 

--  20 


--  10 

1   -  20 

--  10 

1   --  20 


NV  number  of  selected  frequencies   integer 

NR  number  of  fault  identification 
regions 


lower  B.P.V.  of  x-|  for  fault 
identification  region  R-l 

upper  B.P.V.  of  x-j  for  fault 
identification  region  R-l 


lower 
upper 


x2 
x2 


lower 
upper 


nv 


nv 


lower     "    x-|    " 

identification  region  R-2 

upper     "    X] 

identification  region  R-2 


decimal 


2NV+1 


I  -  10 

II  -  20 


lower 


upper 


nv 


nv 
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2NV+2 


1       —     10  lower  "        x'-| 

identification  region  R-3 


11 


20 


upper 


identification  region  R-3 


NR-NV+1 

NR-NV+2 

1 

— 

10 

NM     number  of  measurements 

integer- 

NR-NV+3 

1 

-- 

10 

magnitude  value 

of  X-| 

decimal 

11 

-- 

20 

ii 

x2 

ii 

21 

NV 

30 

ii 
1. 

x3 

ii 

NR-NV+4 

1 

— 

10 

magnitude  value 

Of   X] 

ii 

11 

— 

20 

ii 

*2 

H 

NV 


1. 


NR-NV+2+NM 
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PROGRAM  2 

Card 

Ca 

rd 

number 

col 

umn 

1 

1 

— 

10 

11 

-- 

20 

2 

1 

-- 

10 

11 

-- 

20 

21 

__ 

30 

Information 

NI     number  of  M  rows 
NR     number  of  M  columns 
entry  of  col-1,  row-1   of  M 
row-2       " 
row- 3       " 


Form 


integer 


decimal 


row-m 
1       --     10  entry  of  col -2,  row-1   of  M 

11-20  "  row-2       " 

row-ni 


NR+1 

• 

NR+2 

1 

— 

10 

NM     number  of  measurements 

integer 

NR+3 

1 

— 

10 

magnitude  value  of  x-| 

decimal 

11 

NV 

20 

x2 
1. 

ii 

NR+4 

1 

— 

10 

magnitude  value  of  x-| 

n 

11 

— 

20 

'                 *2 

M 

NV 


1. 


NR+2+NM 


*  Open  end  value  may  be  used  +100.   for  upper  boundary  and  0.   for 
lower  boundary. 

*  Appropriate  modification  may  be  done. 
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